Hall effect in strongly correlated low dimensional systems 
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We investigate the Hall effect in a quasi one-dimensional system made of weakly coupled Luttinger 
liquids at half filling. Using a memory function approach, we compute the Hall coefficient as a 
function of temperature and frequency in the presence of umklapp scattering. We find a power-law 
correction to the free-fermion value (band value), with an exponent depending on the Luttinger 
parameter K p . At sufficiently high temperature or frequency the Hall coefficient approaches the 
band value. 

PACS numbers: 71.10.Pm, 71.10.Hf, 72.80.-r, 72.80.Le 
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I. INTRODUCTION 

The Hall effect has been continuously playing an im- 
portant role in experimental condensed-matter research, 
mostly because the interpretation of Hall measurements 
is rather simple in classical Fermi systems.— In such ma- 
terials the Hall coefficient is a remarkably robust prop- 
erty, which is unaffected by interactions and only depends 
upon the shape of the Fermi surface and the sign of the 
charge carriers. Deviations from this simple behavior are 
generally taken as evidence for the onset of strong corre- 
lations and a failure of the Fermi-liquid (FL) paradigm. 2 
Several authors have investigated the Hall effect in three- 
and two-dimensional FL^^ but the question of the role 
of correlations in the Hall effect for low-dimensional sys- 
tems remains largely unexplored. 

In most three-dimensional systems the interactions 
play a secondary role and the FL picture is appropriated 
However, the prominence of interactions increases as the 
dimensionality of the systems decreases and the FL the- 
ory is believed to break down for many two-dimensional 
systems like, e.g., the high-T c cuprate superconductors Z 
In one-dimensional (ID) systems interactions are dom- 
inant, and the FL description must be replaced by the 
Luttinger liquid (LL) theory.-' 9 This theory predicts a 
rich variety of physical phenomena, such as spin-charge 
separation or non-universal temperature dependence of 
the transport properties, 10 many of which have been ob- 
served experimentally. Therefore large deviations from 
the classical Hall effect are expected to occur in quasi- 
one dimensional systems. 

Among the various experimental realizations of low- 
dimensional systems (organic conductors;** carbon 
nanotubeSfi 2 " ultra cold atomic gases^ etc.) the or- 
ganic conductors are good realizations of quasi-f D mate- 
rials. Studies of the longitudinal transport have success- 
fully revealed signatures of LL properties ! 11 : 14 ' 15 Trans- 
port transverse to the chains has given access to the 
dimensional crossover between a pure ID behavior and 
a more conventional high-dimensional onei 10 : 11 ' 16 : 17 To 
probe further the consequences of correlations in these 
compounds, several groups have undertaken the challeng- 
ing measurement of the Hall coefficient i?H(r , )' 18ll9 ' 2 ° l21 
The results, different depending on the direction of the 



applied magnetic field, proved difficult to interpret due to 
a lack of theoretical understanding of this problem. This 
prompted for a detailed theoretical analysis of the Hall 
effect in quasi- ID systems. A first move in this direc- 
tion was reported in Ref. |22| where the Hall coefficient of 
dissipationless weakly-coupled ID interacting chains was 
computed and found to be T-independent and equal to 
the band value. This surprising result shows that in this 
case i?H, unlike other transport properties, is insensitive 
to interactions. However the assumption of dissipation- 
less chains is clearly too crude to be compared with re- 
alistic systems for which a finite resistivity is induced by 
the umklapp interactions^ 

In this work we examine the effect of umklapp scatter- 
ing on the T-dependence of the Hall coefficient in quasi- 
1D conductors. We consider 1/2-fiiled ID chains and 
compute Rb_(T) to leading order in the umklapp scat- 
tering using the memory function approach ] 24 : 25 We find 
that the umklapp processes induce a T-dependent cor- 
rection to the free-fermions value, and this correction de- 
creases with increasing temperature as a power-law with 
an exponent depending on interactions (Fig. [5]). We dis- 
cuss the implications for quasi-lD compounds. 



II. MODEL AND METHOD 

Our model is sketched in Fig.[T] We consider ID chains 
coupled by a hopping amplitude t± supposedly small 
compared to the in-chain kinetic energy. The usual LL 




FIG. 1: Schematics of the model. The chains and the current 
I go along the x-axis, the magnetic field H is applied along 
the 2-axis, and the Hall voltage is measured along the y-axis. 
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model of the ID chains assumes that the electrons have a 
linear dispersion with a velocity vp. For a strictly linear 
band, however, the Hall coefficient vanishes identically 
owing to particlc-holc symmetry A band curvature close 
to the Fermi momenta ±fcp is thus necessary to get a 
finite i?H- We therefore take for the ID chains of Fig. [1] 
the dispersion 



£±(fc) = ±v F (k =p k F ) + a(k =p fc F ) 2 



(1) 



The upper (lower) sign corresponds to right (left) moving 
electrons. Eq. (fTJ) can be regarded as the minimal model 
which gives rise to a Hall effect, while retaining most 
of the formal simplicity of the original LL theory, and 
its wide domain of validity. In particular, this model is 
clearly sufficient at low temperatures (compared to the 
electron bandwidth) since then only electrons close to the 
Fermi points contribute to the conductivities. 

Our purpose is to treat the umklapp term perturba- 
tively. We express the Hamiltonian as Ho + H u where 
H u is the umklapp scattering term and Ho reads 



Ho — 



J jo 

+52 ^jaR^joR^joL^joL 



h.c. 



(2) 



In Eq. ([2]) j is the chain index, T3 is a Pauli matrix, 
and Ajji = JJ A • dl. We choose the Landau gauge 
A y = Hx, such that A/j+i = Hxa y with a y the inter- 
chain spacing. %)A — (ip^ ip L ) is a two-component vector 
composed of right- and left-moving electrons. The sec- 
ond term in Eq. ^ is the band curvature, the third term 
is the forward scattering and the last term corresponds 
to the coupling between the chains. In Eq. ||3J| we have 
omitted the backscattering terms (gi processes) which 
are, for spin rotationally invariant systems, marginally 
irrelevant. 9 We therefore take gi± = g-^n = 0. At 1/2 
filling the umklapp term reads 

Wu = Y J dX J2 {fjoR^j-oR^joL^^L + h - C 
jo 

(3) 

It corresponds to a process where two electrons with op- 
posite spins change direction by absorbing a momentum 
4fcp from the lattice. 

The Hall resistivity p yx relates to the conductivity ten- 
sor through 



Pyx 



(4) 



Here we calculate p yx using a memory function 
approach^ One rewrites the conductivity tensor in 
terms of a memory matrix M(u>) as 



a T {u) = i{lo1 + X (0) [O + iM(uj)} ^(O)} 1 x (0) 



(5) 



where <x T denotes the transpose of <x. The advan- 
tage provided by the memory function is the possibil- 
ity to make finite-order perturbation expansions which 
are singular for the conductivities due to their resonance 
structured This formalism is especially useful for LL 
(Ref. I23I ) and was also used to estimate the Hall coeffi- 
cient in the 2D Hubbard model. '° x(0) is a diagonal ma- 
trix composed of the diamagnetic susceptibilities in each 

direction, x (0) = ( x * (0) J (0) ), with 



dx 



j 



8 2 H 
dAl{x) 



(6) 



A cl =0 



The thermodynamic average (■ • • )o is taken with respect 
to Ho and H is the Hamiltonian of Eq. ([2]) in the presence 
of electric and magnetic fields, A = A cl + A mag . The 
frequency matrix Jl in Eq. ([5]) is defined in terms of the 
equal-time current-current correlators as 5 



a. 



1 



(7) 



From Eq. {5} one can directly express the memory matrix 
M in terms of the conductivity tensor. In the following 
we will only need the off-diagonal term M xy given by 



iM xy (uj) 



<T xx (w)<7yy(u>) + ol y {u) 



Q X y. (8) 



It is then straightforward to rewrite the Hall coefficient 
Rh — Pyx 1 'H as 



Rh(co) = - 



1 



, lim 

*Xi/(0) h^o 



ftxy + iM xy {uj) 
H 



(9) 



III. RESULTS 



From Eqs ^ and ^ we obtain the longitudinal and 
transverse diamagnetic terms as 



X*(0) 



2e 2 v F 



(10a) 



Xy (0) = -2e 2 t x a 2 y j dx(^(x)^(x)e~ ieHa « x + h.c.) 

(10b) 

For evaluating the frequency matrix we write down the 
current operators: 

J x = e I dx ip] a (x) [vfT3 + 2a(-id x )Ti] ip itr (x) 

J jo 

(11a) 

J y = -iet±ay dxJ2 (^j^ j+ i,cr e ' ieA " +1 - h - c - 

jo 

(lib) 
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The expression resulting from Eq. {7J for the frequency 
matrix is then 



2ireatj_alH 
\ — 

Up 



dx (tpQ* (aOV'it (x)t 



-h.C.). 

(12) 

At this stage we can already evaluate the high-frequency 
limit of i?Hj because the memory matrix vanishes as 1/ui 
(Refs H and [H) and thus M xy drops from Eq. © if 
lo — > oo. The effects of the umklapp disappear at high 
frequency, and in this limit one recovers from Eqs (f9T[T2"|) 
the result obtained for dissipationless chains^ namely 
that the Hall coefficient equals the band value i?^: 



o Traa y 



(13) 



In the definition of the memory matrix, Eq. ([8]) , we can 
ignore the terms of order H 2 which do not contribute to 
i?H in Eq. (|9|) . Furthermore we express the conductivities 
in terms of current susceptibilities as a pv — ^\Xp,v — 
which leads to 



<W„(o)] 



iM xy (uj) 



[Xx(0) - Xxx(w)] [Xy(0) ~ Xyy(u)} 



We rewrite this formula at intermediate frequencies, such 
that |x A , At (w)/x/i(0)| is small. In this expansion we use 
the equation of motion of the susceptibilititespl as well 



as the relation [Hq , J M 



-Sl vp J v . For [i = x the latter 



equation is exactly satisfied in our model, while for fJ, = y 
it is only verified in the continuum limit a y — > 0. The 
resulting expression of the memory matrix is 



iM xy (u) 



1 (K x ;K y )-(K x ;K y ) 



X»(0) 



(15) 



where K„ are the residual forces operators defined as 
the part of the Hamiltonian which in the absence of 
magnetic field does not commute with the currents, i.e. 
K p = [W u , Jfj], and (K x ; K y ) stands for the retarded cor- 
relation function of the operators K p . The terms omitted 
in Eq. (fTS)) are either of second order in |x^(w)/x M (0)|, 
of second order in H, or vanish in the continuum limit 
a y -c 0. Using Eqs © and JTT]) we find 



K, 



h.C. 



2ev-pg3 / dx 

■* j° 

iet ± g 3 a y dx^ ^ 

J ja b=L,R 

(^}^}_ <7b 't(}j,-a-,-b^j+l,<T,-b- 
^j-l^i-ob^j-o-b^ja-b) + h.C. 



(16a) 



(16b) 



Note that each of the .ftT's is of first order in 53, hence M xy 
is of order g|. The quantity (K x ;K y ) entering Eq. (Tl£ 




max(T*, A) 



FIG. 2: Correction of the high-temperature/high- frequency 
Hall coefficient Rh by the umklapp scattering in weakly- 
coupled Luttinger liquids. i?H is the value of the Hall co- 
efficient in the absence of umklapp scattering, Eq. (|13[) . and 
W is the electron bandwidth. Our approach breaks down be- 
low some crossover scale (dashed line, see text). In this figure 
we have assumed that A in Eq. (|19p is negative. 



is the real-frequency long- wavelength limit of the corre- 
lator, which we evaluate as 



/3 



(K x ;K y ) = - / dTe niT (T T K x (T)K y {0)} 



(17) 



It is easy to prove that (K x ;K y ) vanishes for H — or 
a = 0, by applying spatial inversion and particle-hole 
symmetry, respectively. Retaining only leading-order 
terms in t± and a, the first nonvanishing contribution 
in Eq. (fTT|) is of order at^g^H, and involves three spa- 
tial and three time integrations, which we were not able 
to perform analytically. Based on a scaling analysis, we 
can nevertheless extract the temperature (or frequency) 
dependence of this contribution (see Appendix [A]) , which 
yields: 



1 



(K x ;K y 



ixAO)Xy(O) 



a g\ max(w, T) 



3K„ 



(18) 



where K p is the LL parameter in the charge sector. In the 
absence of interactions we have K p = 1, while K p < 1 
(K p > 1) for repulsive (attractive) interactions. If the 
interactions are strong and repulsive (K p <C 1) the ex- 
ponent in Eq. (fT8)) changes due to the contraction^ of 
the operators in K x and K yi which gives the relevant 
power-law in this case. 

Together with Eqs (dHJ) and ©, Eq. flU) leads to our 
final expression for the Hall coefficient: 



\7TVf 



3K„ 



(19) 



with W the electron bandwidth. Eq. (|19p shows that 
in 1/2-filled quasi-lD systems the umklapp scattering 



4 



fex, iv\ k\ + 5k, iv\ 
— O * 



it 



.93 



ii 



lit 




k\ + q, ivi + iv 2 — jV 3 + *fi 



FIG. 3: Example of a diagram appearing in Eq. (I17|l at first 
order in t± and for 92 = 0. The full (dashed) lines correspond 
to free right (left) moving fermions, j is the chain index, and 
the arrows represent up and down spins. The magnetic field 
increases the momentum of the electron by Sk — eHa y . 



changes the absolute value of the Hall coefficient with 
respect to the band value, which is only recovered at 
high temperature or frequency. Note that Eq. (|T9|) also 
describes the frequency dependence of i?n provided T is 
replaced by lu. The backscattering term g\ (neglected 
here) could possibly give rise to multiplicative logarith- 
mic corrections to the power law in Eq. (fTS)) . 9 The sign 
of the dimensionless prefactor A can only be determined 
through a complete evaluation of (K x ;K y ) in Eq. (fT7|) . 
and is for the time being unknown. The available exper- 
imental data are consistent with Eq. (fT9|) if one assumes 
that A is negative (see below) , as illustrated in Fig. [5] 

Eq. p9[) would imply that in the non-interacting limit 
K p — > 1 (g-2 — ► 0) the correction to the Hall coefficient be- 
haves as \og(T/W). In order to check this prediction we 
have evaluated the correlator in Eq. (fT7|) for 92 = 0. The 
corresponding diagram sketched in Fig. [3] involves three 
frequency-momentum integrations, which in this case can 
be fully worked out analytically (see Appendix [BJ . The 
resulting expression of i?n at zero frequency is (T W) 



8 \ttvf J & \W 



(20) 



consistently with Eq. (|T9|) . For non- interacting electrons, 
though, we see that the relative correction induced by 
the 1 / 2-filling umklapp is positive at T < W. Since all 
properties are analytic in K p , we can also deduce from 
Eqs HU) and ([20]) that A tends to [24(1 - Kp)}- 1 in the 
limit K p — > 1. Note that Eq. ([2"0")) would also apply to 
models in which 92 ~ 33, such as the Hubbard model, 
while Eq. (|T9)) is valid only when 173 <C #2- 



IV. DISCUSSION AND CONCLUSION 

The result of Eq. HU) shows that in 1/2-filled quasi ID 
systems the umklapp processes induce a correction to the 
free-fermion value (band value) of the Hall coefficient i?n, 
which depends on temperature as a power-law with an 



exponent depending on interactions. At high tempera- 
tures or frequencies, i?H approaches the band value as 
shown in Fig. [21 implying that any fitting of experimen- 
tal data must be done with repect to the value of i?H at 
high temperature or frequency. 

To study the range of validity of our result, one must 
consider that at low temperature the quasi- ID systems 
generally enter either an insulating state characterized by 
a Mott gap A, or a coherent two- or three-dimensional 
phase below a temperature T* controlled by t±; 10 in ei- 
ther case our model of weakly-coupled LL is no longer 
valid, as illustrated in Fig. [2] The variations of i?H be- 
low max(T*,A) can be very pronounced, and depend 
strongly on the details of the materials. When the ground 
state is insulating, for instance, Rr(T) is expected to go 
through a minimum and diverge like e A / T as T — > 0, 
reflecting the exponentially small carrier density. Other 
behaviors, such as a change of sign due to the formation 
of an ordered state or nesting in the FL regime^ can also 
occur. The validity of Eq. (fTl?|) is therefore limited to the 
LL domain: max(T* , A) < max(T, w) < W. 

For the case A > T* , we estimate the change of Rn 
with respect to at the crossover scale A, for a sys- 
tem with 173 <C U, where U is the Coulomb repulsion. 
The umklapp-induced Mott gap in 1/2-filled systems is 
giver£ by A/W ~ [9 3 /(^v F )] x with x = [2(1 - Kp)]- 1 . 
We thus find that the largest correction is ~ [33/ '(ttvf)]* 
and has a universal exponent. On the other hand, i?u 
approaches the asymptotic value R^ quite slowly, and 
according to Eq. (| 19[) a correction of ~ [93/ (^vf)} 2 still 
exists at temperatures comparable to the bandwidth. 

The available Hall data in the TM family and in the 
geometry of the present analysi a 18 ' 20 show a weak cor- 
rection to the free fermion value which depends on tem- 
perature. Some attempts to fit this behavior to a power 
law have been reported.— However the analysis was per- 
formed by fitting Rn(T) to a power law starting at zero 
temperature. As explained above, the proper way to 
analyze the Hall effect in such quasi-lD systems is to 
fit the deviations from the band value starting from 
the high temperature limit. It would be interesting to 
check whether a new analysis of the data would pro- 
vide good agreement with our results. However in these 
compounds both 1/4-filling and 1 /2-filling umklapp pro- 
cesses are present. For the longitudinal transport, the 
1/4-filling contribution dominates^ For the Hall effect, 
the analysis in the presence of 1 / 4- filling umklapp is con- 
siderably more involved, but a crude evaluation of the 
scaling properties of the corresponding memory matrix 
gives also a weak power-law correction with an exponent 
2— WKp+(K p +K p ~ 1 )/2, and thus similar effects, regard- 
less of the dominant umklapp. The observed data is thus 
consistent with the expected corrections coming from LL 
behavior. However more work, both experimental and 
theoretical, is needed for the TM family because of this 
additional complication, and to understand the data in a 
different geometry where no temperature dependence is 
observed^ 
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Our result Eq. (|T9|) is however directly relevant for 
1 /2-filled organic conductors such as (TTM-TTP)I 3 
and (DMTSA)BF4.— Hall measurements for these com- 
pounds still remain to be performed. Comparison of 
the Hall effect in these compounds with the one in 
1 / 4- filled non-dimerized system o 11 ^ 27 for which only f /4- 
filling umklapp is present, could also help in understand- 
ing the dominant processes for the TM family. 
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APPENDIX A: SCALING ANALYSIS 



In order to establish Eq. (JT8J) we evaluate the correlator 
(K x ;K y ) to first order in a and t±. Let's denote by TL a 
the curvature [second term in Eq. ((2])] , by TL j_ the inter- 
chain hopping [fourth term in Eq. l[2])]. and by TLid the 
remaining part of the Hamiltonian, Hid = H-o—Ti-a^Ti^. 
Standard perturbation theory yields 



(K x ;K y ) 



dr e 



iVtr 



dr\dT 2 



(T T K X (r )K y (0)H_l (ri)H a (T2)) (AI) 

where the average has to be taken with respect to TiiD- 
The latter corresponds to a ID chain and can be easily 
bosonizedi 9 - 



(A2) 
(A3) 



HlD — "Hp + T~(-a 

U v = {u v K v [V9 v {x)] 2 + |i[V^( a ;)] 2 | 

where v = p(a) denotes the charge (spin) degrees of free- 
dom, u v is a velocity, K v a dimensionless parameter, and 
9 V and 4> v are bosonic fields. In our case we have K p < 1 
and K a = 1. The fields 6 V and <f> v and the fermionic 
fields ip are related by 

Abk-px 



=e ^ 1 



,(x)-8 p (x)+cr[b4, tr (x)-e t ,(x)]} 



(A4) 



with b = +1(— 1) corresponding to right (left) moving 
fermions, and a a cutoff of the order of the in-chain lattice 
spacing. With the help of Eq. (|A4j) we bosonize each 
operator in Eq. (|A1[) : 

4iev F g 3 



K x 



(2™) 



dx sin 



ix ejji ( 



(A5) 



— ieHa 



iet±g 3 a y x - x - 

(27ro) 2 ^ ^ 
V ' (j,f) *b 

e ^,{ib<t, p (x)-c ]]l e p (x)-a[b<t,Jx)+c ]] ,e,(xj\} ] 
^={b^ p (x)+e jj ,e p (x)+a[b4, l7 (x)+e jj ,e^(x)]} jl 



+h.c.) 



(A6) 



where j and j' are neighboring chains and e^y 
(— 1) J ~ J . For the hopping term we have 



Hi 



— Y [ dx(e- ieHa y x 
lira^ J V 



jab 

e ^{H P (x)-e p (x)+a[b4> a (x)-e a (x)]} ] 

e -^{H P (x)-e p (x)+a[b(t> tr (x)~e a (x)]} j+1 

-h.c. 



and for the band curvature term we take^ 
na 2iTa / 2 ■ 



(A8) 



Next we use the identity 9 



n < rn 



( JJ^[A„0(r»)+B n9 (r„)] )=exp | 
n 

- (A n A m K + B n B m K~ x )F x {r n - r m ) 

+ (A n B m + B n A m )F 2 (r n - r m )} (A9) 

where r = (x, ut), the notation means that the sum is 
restricted to those terms for which J2 n ^« = B n = 0, 
and i*x,a & r e universal functions. The resulting expression 
for the correlator in Eq. (|A1|) is 

(K x ;K y ) ~H J d 2 rd 2 r 1 d 2 r 2 e- 3K " Fl{r) \r\ 



e -K p F 1 (r-r 1 ) e ±(K p -K p 1 -2)F 1 (r 1 )_ 



\ r 2\ 



(A10) 



The factor |r| results from the linearization in H, and 
we have discarded all factors involving the F 2 function, 
since they correspond to angular integrals of the r vari- 
ables and therefore do not contribute to the scaling di- 
mension. At distances much larger than the cutoff a we 
have e~ AFl ^ ~ (a/|r|)' 4 , and therefore we find the high 
temperature, high frequency behavior as 

(K x ;K y ) ~ ffmax^T)- 3 ^'"-^-^ 1 ). (All) 

We follow the same procedure for the diamagnetic term 
Xj/(0) — however at zeroth order in a and H — and find 



(AI2) 



Combining these expressions and collecting the relevant 
prefactors we deduce Eq. ([18]) . 



APPENDIX B: HALL COEFFICIENT WITHOUT 
FORWARD SCATTERING 

Here we provide the derivation of Eq. (|2U)l , which gives 
i?H to leading order in g 3 but in the absence of for- 
ward scattering (g 2 = 0). Using Eqs §§§ and (fT3"]) we 
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can express the zero-frequency Hall coefficient in terms 
of R^ and Re M xy (iO + ). We then perform a Kramers- 
Kronig transform, insert the free-fermion values of the 
diamagnetic susceptibilities, Xx{0) — —2e 2 v-p / (ira y ) and 
Xy(0) = — ^t^ciy / (itvf) , and use Eq. (fl~5)) to arrive at 



Rh(0) = K 



8e 3 at^a y H 



Im 



K, 



y/°\in->Lu+io+ 



(Bl) 



where (K x ;K y )o is to be evaluated to first order in H. 



The lo = term in Eq. (|15p disappears due to the princi- 
pal part in the u> integral in Eq. (|Blj) . From Eq. (| 16|) one 
sees that (K x \ Ky)o involves 8 fermion fields and can be 
represented by diagrams like the one displayed in Fig. [3] 



There are 32 different diagrams, but all of them can be 
expressed in terms of only one function A(iQ, H), whose 
expression is given by the diagram in Fig. [3] We thus 
obtain 



i? H (0) = \l 



Hal 



(III! 



Im [A'(u> + i0 + ) - A(-co - i0 + )] \ (B2) 



where A'{iO) = dA(itt, H)/dH\H=o and we have pulled 
all prefactors from Eq. p^|) . as well as a factor t± from 
the diagram, out of the definition of function A. The 
explicit expression of A' is 



A'(iQ) = 



(2 



TP)* / 



dkidk2<lq 



dMh) 1 



I 



1 



1 



iv2 - £,+ (fe) 1V3 - £- (&2 - q) iv\ + iv 2 — iv-i + iQ — ^ (fci + q) 



(B3) 



The frequency summations in Eq. (|B3j) are elementary, 
and the various momentum integrals can also be evalu- 
ated analytically to first order in a, yielding 



Ru(0) = R° H 



I - 



1 

16 \ttvf 
(/3 W /4) 2 



da; 



sinh 2 (/3w/4) 



tanh(/3w/4) sinh 2 (^w/4) 



(B4) 



The remaining energy integral is divergent and must be 
regularized. Cutting the integral at the bandwidth W 
and assuming T <C W we obtain the asymptotic behavior 
given in Eq. (f2"TJ|) . 
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